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INTRODUCTION 


The  Loma  Prieta  Earthquake  of  October  17,  1989  caused  $125  million 
dollars  of  damage  to  U.S.  Navy  facilities.  The  predominant  cause  of  damage  was 
liquefaction  of  cohesionless  waterfront  deposits.  This  is  a  major  continuing  problem 
faced  by  the  Navy  which  because  of  mission  requirements  must  be  situated  at  the 
waterfront  often  on  marginal  soils.  This  report  will  analyze  the  response  of  soils  at 
Naval  Station,  Treasure  Island  documenting  settlements  observed  and  serving  as  a 
benchmark  to  calibrate  soil  settlement  prediction  algorithms.  Of  significance  to  this 
event  was  the  fact  that  80  percent  of  the  loss  of  life  occurred  in  a  heavily 
concentrated  damaged  area  approximately  50  miles  from  the  fault  rupture  zone. 
The  local  site  soil  effects  are  the  most  striking  feature  of  the  Loma  Prieta  event. 
Reference  1  through  20  were  used  as  part  of  this  study. 

The  Earthquake 

The  earthquake  (References  1,2,3)  occurred  when  a  segment  of  the  San 
Andreas  fault  northeast  of  Santa  Cruz,  California  ruptured  over  a  length  of  28  miles 
producing  a  Richter  local  magnitude,  Ml,  of  7.0  and  an  average  surface  wave 
magnitude,  Ms,  of  7.1.  The  epicenter  was  10  miles  northeast  of  Santa  Cruz  and  20 
miles  south  of  San  Jose.  Figure  1  shows  the  fault  rupture  main  shock  and  the 
numerous  after  shocks  in  both  plan  view  and  in  two  section  views.  Figure  2 
illustrates  a  cross-section  of  the  fault  region  showing  the  dipping  fault  plane  and  the 
hypocenter,  Reference  2.  The  initial  rupture  length  was  estimated  to  be  24  miles. 
The  main  rupture  began  at  a  depth  of  1 1  miles  below  the  earth's  surface  and  near 
the  center  of  what  would  be  the  rupture  plane.  Over  the  next  7  to  10  seconds  the 
rupture  spread  approximately  12  miles  to  the  north  and  12  miles  to  the  south  The 
unusual  middle  location  of  the  hypocenter  within  the  rupture  location  contributed  to 
the  unusually  short  duration  of  the  event.  Approximately  8  to  10  seconds  of  strong 
shaking  was  observed  which  is  considerably  less  than  would  be  expected  from  an 
event  of  this  size.  The  rupture  propagated  towards  the  earth's  surface  but  during  the 
main  event  appears  to  have  stopped  at  a  depth  of  3  to  4  miles. 


The  section  of  fault  rupture  for  the  Loma  Prieta  earthquake  is  the  southern  part  ■ — - 

of  the  section  of  the  San  Andreas  which  ruptured  in  the  1906  magnitude  8+  ~ 
earthquake.  This  part  of  the  fault  contains  a  bend  resulting  in  a  localized 
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compression  zone.  Figure  2  illustrates  the  fault  showing  the  right  lateral 
(horizontal)  and  reverse  ship  (vertical)  motion.  The  right  lateral  offset  was  6.2  feet 
and  the  vertical  offset  was  4.3  feet,  Reference  1.  This  combination  was  a  significant 
cause  of  the  unusually  high  vertical  ground  surface  accelerations.  Figures  3  and  4 
show  the  ratio  of  vertical  acceleration  to  horizontal  acceleration.  In  Figure  3  both 
components  are  considered  and  in  Figure  4  the  average  horizontal  component  is 
used.  These  figures  show  that  in  the  region  near  faulting  vertical  components 
exceed  the  traditional  level  of  50  percent  of  horizontal  components.  A  total  of  51 
aftershocks  of  magnitude  3.0  or  larger  occurred  within  the  first  24  hours  after  the 
main  shock.  Figure  5  shows  the  region  affected  by  the  event.  Several  factors 
combined  to  keep  damage  and  loss  of  life  down: 

The  rupture  occurred  in  a  sparsely  populated  location. 

The  duration  of  strong  motion  was  short. 

The  event  hypocenter  was  centrally  located  and  spread  uniformly. 

NAVAL  STATION  TREASURE  ISLAND 
Geology 

Figure  6  shows  the  faulting  around  the  Treasure  Island  site  and  the  Loma 
Prieta  earthquake  rupture  zone  (References  4,5  and  6).  Treasure  Island  is  a  man¬ 
made  island  constructed  in  the  1930s  and  situated  between  San  Francisco  and 
Oakland  and  attached  to  Yerba  Buena  Island  by  a  short  causeway,  Figure  7, 
Reference  7.  Over  29  million  cubic  yards  of  mostly  fine-to-medium  grained  sand 
was  dredged  from  borrow  areas  in  the  San  Francisco  bay  and  used  as  fill  material 
over  the  Yerba  Buena  Shoals  north  of  Yerba  Buena  Island.  The  bottom  of  the 
shoals  area  varied  in  depth  from  -2  feet  to  -  26  feet  mean  lower  low  water.  About  65 
percent  of  the  bottom  sediments  in  this  area  were  composed  of  sand  and  the 
remainder  was  soft  clay.  A  low  mound  of  rock  was  placed  along  the  perimeter  of 
island  to  act  as  a  retaining  dike  for  the  sand  fill,  Figure  8.  The  fill  material  was 
deposited  hydraulically  by  using  a  pipeline,  by  hopper  and  by  clam  shell  dredge. 
Where  the  depth  of  the  shoals  exceeded  -6  feet  a  bed  of  hydraulic  sand  fill  was 
placed.  The  dikes  were  constructed  such  that  each  succeeding  level  was  placed 
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inward  of  the  lower  dike  and  rested  on  previous  levels  of  hydraulic  fill.  Fill  was 
placed  to  a  level  of  + 13  feet.  Photography  of  the  construction  shows  that  the  dike 
was  constructed  in  segments  starting  at  the  southwest  corner.  A  low  weir  was 
installed  to  allow  water  from  the  hydraulic  dredges  and  soft  mud  displaced  by  the  fill 
to  escape.  A  small  dredge  was  used  to  remove  areas  of  entrapped  mud. 

During  construction  a  500  foot  length  of  the  north  end  of  the  east  perimeter  dike 
settled  10  to  14  feet.  This  failure  area  was  stabilized  by  flattening  the  slope  and 
placing  a  layer  of  sand  beyond  the  dike  toe.  The  north  dike  design  was  modified  by 
first  excavating  a  400-foot  trench  20  to  30  feet  deep  and  backfilling  with  coarse  sand 
as  a  foundation  for  the  rock  dike,  Reference  7. 

Soil  Conditions 

The  ground  level  of  Treasure  Island  varies  from  + 10.5  to  + 16  feet  above  mean 
lower  low  water  level  with  a  water  table  between  +6  feet  and  +0  feet.  Water  levels 
are  affected  by  the  permeability  of  the  sand  fill  and  vary  with  proximity  to  the 
perimeter  dike.  The  ground  water  levels  at  the  center  of  the  island  are  less  affected 
by  the  tide  and  vary  from  5  to  8  feet  below  the  ground  surface.  During  the 
earthquake  it  is  thought  the  water  table  near  the  island  perimeter  was  at  +3  feet. 
Numerous  explorations  have  been  made  by  drill  and  sample  or  cone  penetration 
testing,  Figures  9  and  10,  References  4,  5  and  6.  Subsurface  materials  can  be 
divided  into  four  strata: 

1)  Loose  to  medium  dense  hydraulically  placed  sand  fill; 

2)  Loose  to  medium  dense  native  material,  Yerba  Buena  Shoals  sands  and 
medium  stiff  native  clays, 

3)  Recent  Bay  sediments  of  medium  stiff  olive  gray  silty  clay  (Bay  Mud)  but 
containing  some  soft  clay  zones;  and, 

4)  Older  Bay  mud  sediments  consisting  of  brownish  and  greenish  gray  very  stiff 
sandy,  silty  and/or  peaty  clays  and  dense  sands. 
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The  fill  material  is  a  fine  to  medium  grained  sand  and  has  gradations  ranging 
from  well  to  poorly  graded.  It  contains  varying  amounts  of  gravel,  silt  and  clay.  In 
general  the  fill  material  is  of  lower  density,  looser,  has  a  lower  penetration 
resistance  and  lower  shell  content  than  the  native  shoals  sands  and  clays.  However 
from  an  engineering  perspective  the  fill  and  native  materials  are  thought  sufficiently 
the  same  from  an  earthquake  performance  basis  so  as  to  be  classified  as  a  common 
strata.  Thickness  contours  of  the  fill  and  native  shoals  strata  are  shown  in  Figure  1 1 
and  is  seen  to  vary  between  30  to  50  feet  for  the  most  part.  Figure  12  gives  cross 
sections  of  Standard  Penetration  Test  corrected  blowcounts  for  this  strata.  The 
blowcount  data  is  reported  from  recent  data  using  appropriate  penetration 
procedures  (rotary  wash  drilling  techniques,  appropriate  hammer  size  and  drop  etc.) 
and  from  Cone  Penetiometer  Test  probes  from  which  data  was  converted  to 
blowcounts  utilizing  site  specific  correlations.  Corrected  data  from  previous 
investigations  has  also  been  included.  The  blowcounts  show  the  material  to  be  loose 
to  medium  dense  and  susceptible  to  liquefaction  under  seismic  shaking.  Note  the 
extensive  shaded  region  containing  blowcounts  less  than  10.  Values  less  then  10 
would  have  a  high  probability  of  liquefaction  under  moderate  earthquakes  and  are 
high  hazard  regions  regardless  of  what  procedure  of  analysis  is  used. 

Below  the  fill  and  native  material  layer  is  a  layer  of  Bay  Mud  composed  of  a 
medium  plastic  silty  clay  with  interbedded  regions  of  sand  and  silt.  These  recent 
Bay  sediments  vary  in  thickness  between  10  and  120  feet.  The  layers  were  deposited 
in  a  marine  environment  and  are  normally  consolidated. 

The  older  Bay  sediments  of  Pleistocene  age  are  generally  stiff  to  sandy,  silty  and 
or  peaty  clays  that  extend  down  to  the  Franciscan  bedrock.  The  layer  is  lightly  to 
moderately  overconsolidated.  Bedrock  is  at  a  depth  of  about  280  feet  below  ground 
surface.  Thickness  of  the  older  Bay  sediments  is  estimated  to  range  from  20  to  170 
feet. 

Figure  13  shows  a  plan  view  of  Treasure  Island  where  post  earthquake  testing 
was  performed  by  Hryciw  et.  al..  Reference  5.  Based  on  seismic  cone  penetration 
test  shear  wave  propagation  velocities  for  the  fill  materials  and  the  recent  Bay  Muds 
were  determined,  Figures  14  and  15.  Equations  of  best  fit  were  computed.  Shear 
modulus  and  damping  as  functions  of  shear  strain  were  determined  for  use  in 
analysis,  Figure  16. 
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Ground  Motion 


Strong  motion  recordings  were  obtained  from  an  instrument  on  Treasure  Island, 
Figure  17,  References  1,  2  and  5.  The  peak  horizontal  ground  acceleration 
components  from  the  main  shock  were  0.1 6g  and  0.1  Og.  A  significant  factor  in  the 
Loma  Prieta  earthquake  was  the  amplification  of  ground  motion  in  areas  underlain 
be  thick  deposits  of  Bay  sediments.  Treasure  Island  falls  within  this  observation 
especially  in  comparison  with  recordings  on  nearby  Yerba  Buena  where  the  peak 
horizontal  acceleration  recorded  on  a  rock  site  were  about  three  times  less  than 
those  on  Treasure  Island.  This  will  be  discussed  in  the  following  section. 
Observation  of  the  Treasure  Island  record  shows  that  at  about  15  seconds  after  the 
start  of  recording,  the  ground  motion  was  subdued;  this  was  probably  caused  by  the 
occurrence  of  subsurface  liquefaction.  Liquefaction  occurred  after  about  4  or  5 
"cycles"  of  shaking  after  about  5  seconds  of  strong  motion.  Sand  boils  were  observed 
at  numerous  locations  and  bayward  lateral  spreading  occurred  with  associated 
settlements.  Ground  cracking  was  visible  with  individual  cracks  as  wide  as  6  inches. 
Overall  lateral  spreading  of  1  foot  was  estimated.  Ground  survey  measurements 
indicate  that  settlements  of  2  to  6  inches  occurred  variably  across  the  island  and  that 
some  areas  had  as  much  as  10  to  12  inches  of  settlement.  The  liquefaction  related 
deformations  resulted  in  damage  to  several  structures  and  numerous  broken 
underground  utility  lines. 

Yerba  Buena 

Yerba  Buena  Island  is  a  large  rocky  outcrop.  Figure  18  shows  the  accelerogram 
components  recorded  for  the  main  shock.  Note  that  the  horizontal  components 
were  0.068g  and  0.03  lg,  both  significantly  less  than  those  on  Treasure  Island. 
Figures  19  and  20  show  the  5  percent  damped  response  spectra  for  the  two  sites  for 
both  horizontal  components  of  motion  illustrating  the  spectral  amplification.  Note 
the  significant  soft  soil  amplification  shown  in  Figure  19  for  Treasure  Island  and 
compare  it  to  Figure  20,  Reference  2. 
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EARTHQUAKE  DAMAGE 


General  Ground  Deformation 

There  was  extensive  evidence  of  settlement  and  differential  settlement 
adjacent  to  or  inside  buildings.  The  region  around  Buildings  2  and  3  which  was 
constructed  on  piles  settled  on  the  order  of  4  to  6  inches  while  the  structures 
remained  at  their  original  elevation.  There  was  evidence  of  extensive  liquefaction 
by  the  presence  of  sand  boils  up  to  20  feet  in  diameter.  The  large  sand  boils  with 
the  accompanying  settlement  and  expulsion  of  water  created  areas  of  ponded  water. 

Figure  13  shows  locations  of  test  performed  by  Hryciw  (Reference  5).  At  the 
UM1  location  no  damage  was  seen.  Some  liquefaction  occurred  in  adjacent  inland 
areas  near  UM3  and  subsidence  of  the  retaining  dike  was  noted  at  UM12.  During 
construction  the  north  eastern  section  of  the  dike  settled  extensively  and  the  dike 
design  was  modified.  At  this  location,  a  400  foot  long  trench  was  dug  to  a  depth  of 
of  20  to  30  feet  and  was  backfilled  with  sand. 

Lateral  spreading  and  settlement  of  the  dike  occurred  in  a  number  of 
locations.  The  maximum  settlement  at  the  top  of  the  dike  was  nearly  2  feet.  Lateral 
spreading  cracking  was  widely  observed.  It  is  estimated  that  8  inches  of  the 
settlement  was  caused  by  shaking-induced  compaction  and  the  remainder  was 
lateral  spreading.  A  2500  foot  long  crack  opened  near  the  east  side  of  the  island 
and  passed  through  Building  7.  The  floor  slab  cracked  and  a  major  sand  inflow 
occurred  filling  the  ground  floor  of  the  building  with  sand  to  a  depth  of  6  inches  in 
one  area.  The  area  between  3rd  and  5th  Streets  and  near  11th  Street  exhibited 
lateral  spreading  movements  having  cracks  of  0.5  inches.  Summation  of  horizontal 
crack  widths  indicates  a  total  bayward  movement  of  6  to  12  inches  may  have 
occurred  at  the  dike  area.  Survey  data  shows  that  at  Avenue  N  and  3rd  Avenue 
approximately  10  inches  of  lateral  movement  occurred  primarily  eastward.  Most 
spreading  movement  was  restricted  to  an  area  within  a  distance  of  200  feet  from  the 
perimeter  dike.  The  lateral  spreading  occurred  only  during  the  period  of  actual 
earthquake  motion  and  did  not  continue  to  move  afterward. 

Ground  surface  settlements  of  up  to  12  inches  occurred.  Examination  of 
subsurface  conditions  indicates  that  the  liquefied  zones  may  have  extended  to  within 
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about  5  feet  of  the  ground  surface.  Six  buildings  suffered  significant  damage.  Several 
structures  of  the  1930's  vintage  suffered  cracking,  separation  of  floor  slabs  and 
minor  sand  inflows.  More  than  40  underground  utility  breaks  occurred.  Pre-  and 
post-earthquake  survey  data  were  used  to  construct  Figure  21  showing  the  shaking 
induced  compaction  settlements  as  reported  by  Egan  et.  al,  References  4  and  8. 
They  used  blowcount  data  to  estimate  settlements  to  expand  the  limited  settlement 
survey  data.  They  initially  used  a  procedure  by  Tokimatsu  and  Seed  (Reference  14) 
which  was  found  to  over  estimate  settlements  by  a  factor  of  two.  A  modified 
procedure  was  developed  which  incorporated  the  liquefaction  assessment  models  of 
Liao  and  others  (Reference  15)  and  the  influence  of  grain  size  on  post-cyclic 
volumetric  strain  after  Lee  and  Albaisa  (Reference  16).  Figure  21  shows  elevation 
changes  at  survey  monuments  or  reference  points.  There  was  good  agreement 
between  the  estimated  settlements  and  observed  survey  settlement  data.  Figure  21 
shows  that  for  the  thickest  zones  of  liquefied  soil,  settlements  of  8  to  10  inches  may 
have  occurred. 

Most  buildings  on  Treasure  Island  suffered  no  or  minor  damage  which  was 
limited  to  minor  cracking  or  differential  settlement.  Several  buildings  near  the 
perimeter  dike  or  areas  of  significant  lateral  spreading  did  experience  greater 
damage.  These  areas  experienced  6  to  10  inches  of  settlement.  Buildings  which 
were  situated  in  areas  which  settled  less  than  6  inches  generally  experienced  only 
minor  damage.  Figure  22  summarizes  damage. 

Forty-four  breaks  in  utility  pipelines  included  28  fire  and  freshwater  lines  of  steel 
or  asbestos  cement,  10  sewage  lines  of  vitrified  clay  and  6  welded-steel  gas  lines. 
Many  of  the  breaks  occurred  near  the  dike  in  areas  of  high  lateral  spreading.  Crude 
estimates  of  lateral  spreading  required  to  cause  failure  are: 


Type 

Pipe  Diameter 

Spreading  to  Induce 
Failure 

Steel  or  Asbestos  Cement 

0  to  4  in 

1  inch 

Steel  or  Asbestos  Cement 

12  to  16  in 

6  to  12  inches 

vitrified  clay  pipe 

1/4  inch 
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Densified  Soil 


The  area  by  the  wharf  at  the  south-eastern  comer  of  the  island  was  densified  by 
Terraprobe.  It  showed  no  signs  of  liquefaction,  spreading  or  settlement.  In  areas 
adjacent  to  the  improved  area  major  sand  boils  were  observed  so  it  may  reasonably 
be  concluded  that  the  site  improvement  kept  the  area  from  liquefying.  Lateral 
spreading  was  not  observed  adjacent  to  the  north  section  of  the  dike  which  was 
modified  during  construction  and  placed  on  a  foundation  of  sand,  Reference  13.  In 
this  area  the  dike  subsided  more  than  in  adjacent  areas. 

IMPROVED  SITE  RESPONSE 

Medical  Dental  Building 

The  Medical/Dental  building  was  under  construction  at  the  time  of  the  Loma 
Prieta  earthquake.  The  subsurface  section  along  the  east-west  centerline  of  the 
project  site  is  shown  in  Figure  23.  Note  that  the  soil  to  a  depth  of  31  to  43  feet  is 
loose  to  medium  dense  hydraulically  placed  sand  fill,  generally  fine  to  medium 
grained  containing  less  than  10  percent  material  finer  than  the  #200  sieve.  The 
zone  contains  occasional  thin  layers  of  soft  compressible  silt.  Beneath  this  layer  is  a 
soft  to  medium  stiff  clayey  silt  (Bay  Mud)  layer  of  30  foot  depth.  The  Bay  Mud  is 
underlain  by  alternating  layers  of  dense  to  very  dense  sands  and  stiff  to  hard  clays, 
Reference  13. 

The  structure  was  to  consist  of  a  two  story  steel  frame  structure.  Typical  column 
spacing  was  30  feet  and  the  dead  plus  live  column  loads  are  up  to  230  kips.  The  site 
is  potentially  liquefiable.  The  structure  could  have  been  supported  on  piles  or 
spread  footings  with  site  densification.  It  was  decided  to  density  the  upper  layer  of 
sand  fill  to  a  minimum  relative  density  of  75  percent  beneath  the  building  extending 
to  a  distance  20  feet  beyond  the  building.  Vibro-replacement  technique  using  gravel 
backfill  was  chosen  as  the  method  for  densification.  Tests  were  conducted  to 
determine  probe  spacing.  It  was  decided  to  use  10  foot  probe  spacing  to  a  depth  of 
22  feet.  Since  there  were  zones  of  silt  and  clay  in  the  sand  fill,  the  specifications  did 
not  require  improvement  in  areas  where  the  Cone  Penetration  Test  friction  ratio 
was  greater  than  2  percent.  The  uppermost  few  feet  were  to  be  compacted  to  95 
percent  relative  compaction  using  surface  compactors.  A  study  of  pre  and  post 
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densification  penetration  test  showed  that  the  upper  10  feet  of  sand  was  already 
dense  and  was  not  improved.  The  lower  level  of  the  layer  between  depths  of  10  to 
22  feet  which  was  initially  loose  was  densified.  Using  a  volume  calculation  method  it 
was  estimated  that  the  upper  22  feet  were  densified  to  a  relative  density  range 
between  77  and  80.  The  silty  sand  fill  interbedded  with  zones  of  silt  and  clay  which 
underlay  the  upper  zone  to  a  depth  of  40  feet  was  not  adequately  densified;  no 
attempt  was  made  to  density  soils  in  the  22  to  40  foot  range.  At  the  time  of  the 
earthquake  about  40  percent  of  the  footings  were  built  and  two  22-foot  deep 
elevator  shafts  were  excavated. 

The  bottom  8  feet  of  the  22-foot  deep  elevator  shafts  were  filled  with  sand  from 
liquefaction  flows  of  the  untreated  sand  fill  between  a  depth  of  22  to  40  feet.  Total 
settlement  of  the  site  could  not  be  determined  since  the  benchmarks  settled, 
differential  settlement  of  the  footings  of  0.88  inches  over  a  distance  of  180  feet  was 
noted.  It  appears  that  liquefaction  did  not  occur  in  the  upper  22  feet. 

Building  450 

Building  450  is  a  three  story  steel  frame  office  building  built  in  1967.  The 
structure  is  actually  composed  of  two  buildings,  the  first  160  feet  by  160  feet  in  plan 
and  the  other  54  by  124  feet  in  plan  located  100  feet  northwest  of  the  first.  The 
buildings  have  concrete  walls  and  floors;  typical  dead  plus  live  column  loads  are  250 
to  300  kips.  The  subsurface  section  along  the  north-south  project  centerline  is 
shown  in  Figure  24.  The  profile  is  similar  to  that  in  Figure  23  with  the  upper  30  feet 
being  loose  to  medium  dense  hydraulically  placed  sand  fill  underlain  by  a  layer  of 
soft  to  medium  stiff  silty  clay  approximately  20  feet  thick  (Bay  Mud).  Note  the 
sloping  bedrock  formation,  Figure  24.  It  was  decided  to  use  densification  since 
some  piles  would  be  end  bearing  on  bedrock  and  some  piles  would  only  rely  on 
friction.  Conventional  concrete  piles  were  viewed  as  a  high  cost  alternative  lacking 
lateral  stability  if  liquefaction  occurred.  Sand  compaction  piles  spaced  4  feet  on 
centers  beneath  footings  and  5  feet  on  centers  beneath  floor  slabs  was  chosen  as  the 
method  for  site  improvement.  The  specifications  called  for  a  minimum  relative 
density  of  75  percent  in  the  sand  fill  to  a  depth  of  30  feet  beneath  the  footings  and 
65  percent  to  a  depth  of  30  feet  beneath  the  floor  area  and  to  a  distance  of  10  feet 
beyond  the  building  perimeter.  The  mandrel  used  for  the  sand  piles  was  a  14-inch 
diameter  steel  casing  fitted  with  a  loose  steel  bottom  plate.  The  mandrel  was  driven 
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to  the  required  depth  and  backfilled  with  coarse  sand.  The  top  end  of  the  mandrel 
was  closed  and  a  minimum  air  pressure  of  100  psi  was  applied  to  the  column  of  sand 
as  the  mandrel  was  withdrawn.  Figure  25  shows  Standard  Penetration  Test 
blowcounts  before  and  after  densification.  While  the  average  relative  density 
exceeded  minimum  requirements  and  the  average  overall  densification  was 
adequate,  there  were  isolated  zones  of  silt  and  clay  which  were  not  densified  to 
minimum  requirements,  Reference  13. 

There  was  no  evidence  of  damage  to  this  structure.  Some  lateral  spreading,  sand 
boils  and  localized  settlement  were  observed  outside  the  improved  area. 

Buildings  487, 488  and  489 

The  three  story  buildings  with  exterior  and  interior  concrete  block  walls,  precast 
concrete  floor  slabs,  and  concrete  slab  on  grade  first  floors  were  constructed  in 
1973.  Typical  dead  plus  live  loads  for  the  bearing  walls  are  in  the  range  of  5  to  7 
kips  per  foot.  The  building  does  not  have  columns  or  independent  footings.  Figure 
26  shows  the  subsurface  section  along  the  east-west  centerline  of  the  project  site. 
The  upper  24  to  33  feet  of  soil  are  very  loose  to  medium  dense  hydraulically  placed 
sand  fill  with  occasional  layers  of  soft  compressible  silts  and  clays.  Below  this  layer 
is  a  zone  of  dense  sand  underlain  by  a  4  foot  thick  layer  of  soft  silty  clay  (Bay  Mud). 
The  structure  rests  on  spread  footings  bearing  on  densified  sand  fill.  The 
specifications  required  that  the  sand  fill  be  improved  to  a  relative  density  of  75 
percent  to  a  depth  of  30  feet  beneath  the  buildings  to  a  distance  of  10  feet  beyond 
each  buildings  perimeter.  Vibrocompaction  was  used  for  densification.  The 
specification  called  for  the  vibrator  to  be  inserted  at  each  compaction  point  to  a 
depth  of  30  feet  below  grade  and  maintained  at  that  depth  for  a  period  of  one 
minute  and  then  withdrawn  at  a  rate  of  not  more  than  1  foot  per  minute.  Crushed 
rock  was  continuously  placed  around  the  vibrator  during  the  process.  The  upper  1 
foot  of  sand  was  compacted  to  95  percent  relative  compaction  using  a  vibratory 
compactor.  Standard  Penetration  Tests  taken  indicate  that  relative  densities 
equaled  or  exceeded  specification  minimums. 

Building  487  experienced  minor  cracking  in  the  concrete  floor  slab  caused  by 
differential  settlement.  The  other  buildings  did  not  sustain  any  damage. 
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Buildings  452  and  453 


Buildings  452  and  453  consist  of  six  wings  four  stories  high  radiating  from  a 
central  core.  Building  453,  built  in  1969  has  its  core  loads  carried  by  circumferential 
and  radial  walls;  dead  plus  live  loads  are  approximately  7  kips  per  linear  foot  of 
wall.  Figure  27  shows  a  subsurface  profile  along  a  north-south  centerline  of  the 
project  site.  The  upper  45  feet  consists  of  loose  to  medium  dense  hydraulically 
placed  sand  fill  underlain  by  about  20  feet  of  soft  to  medium  stiff  clayey  silt  (Bay 
Mud).  Below  the  Bay  Mud  are  alternating  layers  of  stiff  clays  and  dense  sands. 
Figure  27  shows  blowcounts  and  indicates  that  liquefaction  would  be  expected  as 
evidenced  by  blowcounts  below  10.  The  building  specification  called  for 
densification  to  70  percent  relative  density  to  a  depth  of  30  feet  under  the  building 
extending  10  feet  beyond  the  building  perimeter.  The  upper  4  feet  were  to  be 
excavated,  backfilled  and  compacted  to  a  minimum  of  95  percent  relative 
compaction.  Nonstructural  timber  piles  were  used  as  the  means  for  densification  of 
both  buildings.  The  piles  were  8  inches  in  diameter  at  the  tip  and  12  inches  at  the 
butt,  20  feet  long  and  driven  25  feet  into  the  fill.  The  tops  were  driven  to  a  depth 
below  the  water  table  to  prevent  deterioration.  Unfortunately  data  giving  post- 
driving  density  is  not  available,  Reference  13. 

After  the  earthquake  the  buildings  were  inspected  and  no  major  structural 
damage  was  observed.  At  Building  453  there  was  one  concrete  spall  and  cracks  in 
the  floor  system.  Some  repairs  were  required  for  the  slab  on  grade  which  had 
settled  less  than  3/8  inch.  No  significant  damage  to  Building  452  was  observed. 
Sand  boils  or  foundation  distress  were  not  observed  around  the  building. 

SOFT  SITE  RESPONSE 

Treasure  Island  Site  Parameters 

Hryciw  et.  al.,  Reference  5,  performed  post-earthquake  cone  penetration  tests  to 
determine  shear  wave  velocities.  This  data  was  presented  above.  Figure  14  and  15. 
They  estimate  that  the  shear  wave  velocity  in  the  sand  fill  upper  layer  was  defined 
by  a  best  fit  of  the  data  as: 


Vs  =  150  +  4  z 
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where  z  is  depth  in  meters  and  Vs  is  in  meters  per  second. 

They  estimate  the  shear  wave  velocity  for  the  Bay  Mud  layer  to  be: 

Vs  =  30  z  0-55 

For  the  alternating  layers  of  stiff  clay  and  dense  sands  the  following  is  estimated: 

Vs  =  250  +  2  z 

For  the  purposes  of  constructing  a  site  model  the  following  additional 
information  can  be  used. 

Soil  Depth  Shear  Wave  Velocity 

Silty  sand  100  to  141  ft.  1100  ft./sec. 

Stiff  to  hard  clay  141ft.  1100ft./sec. 

Stiff  to  hard  clay  285  ft.  1400  ft./sec. 

Hryciw  et.  al..  Reference  5,  used  the  computer  program  SHAKE  to  analyze 
the  soil  profiles  at  Treasure  Island.  SHAKE  is  a  frequency  domain  analysis  using 
strain  dependent  elastic  material  parameters.  They  ran  a  parameter  analysis  to 
consider  the  uncertainty  in  specification  of  the  older  bay  deposits  using  the 
measured  rock  accelerogram  recorded  on  Yerba  Buena  as  the  basis  for  bedrock 
motion.  Their  results  which  would  be  classed  as  typical  response  of  this  category  of 
analysis  are  shown  in  Figure  28.  They  compute  a  peak  surface  acceleration  of  0.18  g 
in  comparison  with  the  recorded  value  of  0.16  g.  The  recorded  spectra  exceeds  the 
computed  spectra  which  is  explained  by  Hryciw  that  the  SHAKE  analysis  does  not 
take  into  account  the  softening  of  the  upper  layers  with  the  onset  of  liquefaction. 
Their  analysis  shows  that  the  upper  layer  most  amplified  the  motion.  Figure  29. 
They  correlated  increase  in  computed  surface  acceleration  with  regions  of  observed 
increase  in  seismic  effects. 

Seed  et.  al.,  Reference  2,  analyze  the  same  site  using  the  following: 
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Soil  Type 
Loose  Sandy  Fill 
Loose  silty  sand 
Silty  clay  Bay  Mud 
Dense  sand  &  silty  sand 
Stiff  to  hard  clay 


Depth 
0  to  30  ft. 

30  to  45  ft. 

45  to  100  ft. 
100  to  140  ft. 
285  ft. 


Shear  Wave  Velocity 
500  to  600  ft./sec. 
550  to  650  ft./sec. 
500  to  700  ft./sec. 
1100  ft./sec. 

1100  to  1400  ft./sec. 


The  calculated  peak  ground  acceleration  of  0.18  g  agrees  with  the  measured  value 
of  0.16  g.  The  results  show  a  spectral  amplification  of  4  to  5  although  the  computed 
spectral  content  is  lower  than  observed,  Figure  30.  Seed  et.  al.  explains  the 
underestimation  of  the  computed  spectra  as  caused  by  surface  waves  generated  by 
the  dipping  of  the  Yerba  Buena  rock  outcrop  beneath  the  alluvium  and  fill  of 
Treasure  Island.  The  late  occurrence  of  liquefaction  would  not  explain  the 
underestimation  of  long  period  spectral  response  (  >  0.15  sec).  They  note  that  peak 
spectral  response  occurs  at  a  period  of  0.6  seconds  and  has  a  secondary  response  at 
1.3  seconds.  The  0.6  seconds  does  not  represent  the  predominant  period  of  this 
deep  soft  site  but  rather  would  be  at  the  1.3  seconds.  The  0.6  second  peak  results 
from  the  site  being  strongly  excited  by  the  input  rock  motion  having  an  energy 
concentration  at  this  period. 

One  should  conclude  that  the  soil  column  analysis  using  SHAKE  is  an  effective 
tool  for  predicting  approximate  levels  of  response.  It  is  a  simple  expedient  tool  to 
use  requiring  data  that  is  usually  easily  available.  It  can  reliably  indicate 
amplification  when  the  site  is  properly  modeled.  However,  it  does  not  account  for 
softening  of  material  properties  or  sloping  bedrock  conditions.  When  applied  to  the 
Treasure  Island  site  it  reliably  warns  the  engineer  that  substantial  amplification  will 
occur.  This  waining  can  be  used  as  a  screening  indicating  additional  analysis  may 
be  needed  when  warranted  by  the  construction  project.  Caution  should  be  exercised 
in  using  code  provisions  which  may  not  be  representative  of  the  worst  case 
conditions  at  a  soft  site.  It  would  be  expected  that  at  high  levels  of  ground  shaking 
significant  nonlinear  response  would  occur  which  would  not  exhibit  the  same  levels 
of  amplification  as  at  low  levels  of  shaking;  however  the  absolute  values  of  motion 
could  be  higher.  Seed  et.  al.  (1992)  shows  results  for  3  additional  sites  where 
SHAKE  was  in  reasonable  agreement  with  observed  spectra.  Figure  31.  They 
conclude  that  nonlinear  methods  in  theory  may  equal  or  exceed  the  accuracy  of 
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linear  methods;  but,  in  practice  the  analyses  often  under-predict  results  because 
overdamping  or  poor  modeling  negate  any  advantage. 


LIQUEFACTION  PREDICTION 
Current  Capabilities 

A  question  of  importance  to  the  Navy  is  whether  we  can  predict  with  high 
reliability  the  locations  on  a  Navy  base  which  will  experience  liquefaction.  In  the 
case  of  Treasure  Island  Standard  Penetration  Blowcounts  of  less  than  10  exist  in 
widespread  areas.  We  have  had  tools  for  over  20  years  which  tell  us  that  areas 
where  blowcounts  are  less  than  10  will  liquefy  during  modest  levels  of  ground 
shaking.  In  1985  Seed  et.  al.  updated  data  using  the  Standard  Penetration  Test 
The  Standard  Penetration  Test  was  standardized  to  an  energy  level  of  60  percent  of 
the  free  fall  energy  of  the  hammer. 

Figure  32  shows  Seed's,  Reference  21,  liquefaction  assessment  chart  where  the 
soil  conditions  are  defined  by  the  Standard  Penetration  Resistance,  (N^)^  and  the 
earthquake  loading  by  average  shear  stress  divided  by  the  effective  vertical 
confining  pressure.  Figure  32  is  based  on  a  magnitude  7.5  earthquake  and  takes 
into  account  the  percentage  of  fine  material  passing  the  number  200  sieve.  Figure 
33  was  generated  to  show  the  effects  of  other  magnitudes.  Both  Figures  32  and  33 
apply  to  level  ground  conditions  free  of  initial  static  shear  stress.  A  soil  element  in 
and  embankment  or  beneath  a  structure  has  an  initial  static  shear  stress  and  has  an 
initial  shear  stress  ratio  of  alpha  defined  as  the  shear  stress  on  the  horizontal  plane 
divided  by  the  effective  vertical  confining  pressure.  A  correction  factor  to  be 
applied  to  the  blowcount  resistance  to  account  for  initial  static  shear  has  been 
developed  by  Hynes  et.  al.  (1988)  based  on  laboratory  tests  on  Folsom  gravels. 
Figure  34.  It  is  assumed  that  the  initial  static  shear  increases  resistance  to 
liquefaction.  Finn  (1991)  points  out  that  "this  increase  applies  to  resistance  to  cyclic 
mobility  rather  than  to  liquefaction,  that  is  to  non-contractive  materials."  Dense 
cohesionless  soils  would  be  expected  to  have  greater  resistance  to  deformation  and 
retain  their  strengths  when  initial  static  shears  are  present.  However  Vaid  and 
Chem,  Reference  19,  show  for  contractive  materials  the  initial  static  shear  decreases 
the  resistance  to  cyclic  loading.  Loose  soils  would  tend  to  experience  an  increase  in 
contractiveness  and  a  strength  drop  after  reaching  peak  strength.  The  effect  of  a 
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soils  looseness  in  terms  of  relative  density  is  shown  in  Figure  35.  When  a  structure 
rests  on  loose  sand  factors  of  much  less  than  1.0  must  be  applied  to  the  shear  stress 
ratio.  This  in  effect  reduces  the  strength  of  the  soil  which  can  be  counted  upon  for 
very  loose  deposits.  This  is  of  critical  significance  to  Navy  waterfront  facilities. 
Increasing  vertical  confining  pressures  also  has  the  effect  of  reducing  the  resistance 
to  cyclic  loading  for  the  very  loose  sands,  Reference  20.  Figure  32  was  developed 
for  clean  sands;  to  apply  Figure  32  to  silty  sands  the  following  corrections  were 
proposed  by  Seed  et.  al..  Reference  18,  to  be  added  to  the  penetration  resistance: 

Fine  Content  Value  to  add 


(blows/ft) 
10%  1 

25%  2 

50%  4 

75%  5 


In  addition  to  the  Standard  penetration  Test,  the  Cone  Penetration  Test  has  been 
related  to  liquefaction.  Reference  22. 

Of  significant  interest  is  that  given  the  occurrence  of  low  blowcounts  what  will  be 
the  resulting  deformation  state  in  a  given  earthquake?  This  is  best  explained  in 
terms  of  the  residual  steady  state  strength  of  the  soil.  Soils  at  a  density  looser  than 
critical  state  are  contractive,  generate  pore  pressures  and  suffer  near  complete  loss 
of  residual  strength,  line  1,  Figure  36.  Cyclic  mobility  is  a  term  applied  to 
cohesionless  soils  that  are  at  a  density  greater  than  that  at  the  critical  void  ratio; 
such  soils  are  expansive  or  less  contractive  than  loose  soils  under  cyclic  loading  and 
exhibit  a  stress-strain  behavior  with  limited  deformation,  line  2  Figure  36.  Under 
loading  dense  soils  tend  to  expand  approaching  the  critical  void  ratio  and  loose  soils 
tend  to  contract  approaching  the  critical  void  ratio.  The  undrained  steady  state 
strength  is  maintained  over  a  large  range  of  strain.  In  1988  Seed  et.  al..  Reference 
18,  developed  an  estimate  of  residual  strength  as  a  function  corrected  blowcounts. 
It  must  be  emphasized  that  this  was  developed  from  case  histories  and  only  is 
intended  to  give  approximate  values. 
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DETAILED  ANALYSIS  OF  SITE 
Site  Analysis 

It  was  of  interest  to  examine  the  Treasure  Island  site  to  explore  reasons  for  the 
site  amplification.  The  data  reported  by  Hryciw  et.  al.  (Reference  5)  was  used  as 
the  starting  point  for  this  investigation.  Additional  information  on  the  site  was 
obtained  from  Reference  23  and  is  shown  in  Figures  37  and  38.  Figure  37  shows  the 
blowcounts  from  a  Standard  Penetration  Test  which  are  noted  to  be  less  than  10  for 
a  significant  portion  of  the  profile.  Figure  38  shows  the  undrained  shear  strength. 

SHAKE  Analysis 

Using  the  same  soil  profile  defined  in  Reference  5,  the  soil  column  was  analyzed 
for  the  Loma  Prieta  earthquake  using  the  SHAKE85  computer  program.  The 
Yerba  Buena  Island  record  was  used  as  the  rock  input  motion  and  average  cyclic 
shear  strain  was  taken  as  0.65  times  the  maximum  shear  strain.  Professor  Rollins, 
Reference  5,  used  SHAKE90  a  more  recent  version  of  the  same  program. 
SHAKE90  uses  strain  dependent  properties  based  on  the  ratio  of  shear  modulus  as 
a  function  of  strain  to  the  maximum  shear  modulus  which  occurs  at  low  strain; 
multiple  functions  can  be  used  for  representation  of  sands  and  clays.  SHAKE85 
uses  a  function  which  is  defined  by  equation;  the  user  is  limited  to  only  1  function 
for  clay.  These  differences  caused  slight  differences  in  the  computed  results.  Figure 

39  shows  the  NCEL  data  which  agrees  very  closely  with  that  of  Reference  5.  This 
establishes  a  benchmark  control  point  from  which  a  detailed  analysis  can  be 
undertaken. 

One  point  of  considerable  interest  is  the  strain  dependent  properties  for  the  Bay 
Mud.  Rollins  in  Reference  24  uses  data  by  Lodde,  Reference  25,  to  define  the 
strain  dependent  shear  modulus  ratio.  Figure  40  shows  a  plot  of  the  Bay  Mud  curve 
compared  with  the  more  normal  values  based  on  data  provided  in  Reference  26.  It 
can  be  seen  that  the  Bay  Mud  has  a  significantly  stiffer  modulus  with  strain.  Figure 

40  also  contains  data  from  Mexico  City,  Reference  27  which  is  very  similar  to  the 
Bay  Mud  behavior.  The  Mexico  City  clays  were  noted  to  be  rather  stiff  at  low  strain. 
Note  that  distant  earthquakes  are  low  strain  events.  The  SHAKE  analysis  was 
repeated  using  the  less  stiff  values  normally  associated  with  clays  as  a  substitute  for 
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the  Bay  Muds.  Figure  41  shows  that  amplification  does  not  occur.  Strains  in  the 
analysis  using  the  Bay  Mud  properties  are  in  the  range  of  0.03  to  0.08  percent  in  the 
Bay  Mud  layers;  this  results  in  an  effective  shear  modulus  of  about  60  percent  of 
maximum  with  damping  in  the  range  of  0.06  to  0.12  of  critical.  However  when 
typical  clay  data  is  used  the  shear  modulus  drops  to  about  10  percent  of  maximum 
and  damping  increases  to  0.08  to  0.15  of  critical.  This  explains  why  the  stiffer  Bay 
Mud  properties  do  not  attenuate  the  motion  as  does  typical  clay.  Reference  24 
evaluates  the  significance  of  shear  wave  velocity  over  a  wide  range  and  concludes 
that  this  is  not  a  factor  of  high  sensitivity.  They  also  conclude  that: 

"the  older  bay  sediments  contributed  very  little  to  the  ground  amplification,  the 
Bay  Muds  contributed  somewhat,  but  by  far  the  greatest  contribution  came  from  the 
fill  material.  However,  this  observation  should  not  be  understood  to  mean  that  the 
fill  sand  is  inherently  more  prone  to  amplification  than  the  Bay  Mud,  but  rather  that 
the  fill  is  under  lower  confining  pressure  and  by  its  surcharging  effect,  provides  the 
Bay  Mud  with  higher  shear  stiffness." 

This  author  disagrees  at  least  in  part  with  this  conclusion.  It  would  appear  from 
looking  at  Figure  39  that  the  amplification  is  occurring  in  the  upper  sand  layer; 
however  as  shown  above  it  is  the  stiff  Bay  Mud  characterization  which  was  shown  to 
cause  amplification.  It  is  true  that  the  confinement  of  the  Bay  Mud  by  the  upper 
sand  layer  contributes  to  the  stiffness  of  the  Bay  Mud.  However  the  Bay  Mud  as 
characterized  by  Lodde  in  Reference  25  shows  it  to  be  stiffer  with  straining  than 
would  be  expected.  This  is  thought  to  be  very  significant.  The  San  Francisco  site 
and  the  Mexico  City  site  both  have  clays  that  are  substantially  stiffer  than  would  be 
expected.  Reference  30  shows  that  the  Plasticity  Index  for  Bay  Muds  is  in  the  range 
of  20  to  40  between  38  and  75  feet.  The  Plasticity  Index  for  Mexico  City  clays  was 
30.  Reference  31  shows  data  documenting  that  the  shear  moduls  is  stiffer  with 
shear  strain  as  the  Plasticity  Index  increases.  This  little  publicized  data  indicates 
that  the  stiffness  of  clay  under  cyclic  loading  should  be  increased  to  account  for  the 
Plasticity  Ratio.  The  Plasticity  Index  is  based  on  the  amount  of  water  required  to 
transform  a  remolded  soil  from  semisolid  to  a  liquid  state.  It  is  a  function  only  of 
the  size  shape  and  mineralogy  of  the  soil  particles  and  the  pore  water. 

The  effect  of  loose  sand  deposits  is  also  clearly  of  interest  since  it  is  believed 
that  these  can  contribute  to  amplification.  This  will  be  explored  in  the  next  section. 
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Effects  of  Loose  Surface  Deposits 


To  study  the  problem  of  whether  loose  sand  deposits  contribute  to  amplification 
a  generic  soil  profile  was  constructed  composed  of  layers  of  sand.  The  upper  70  feet 
was  assumed  to  be  composed  of  sand.  The  Pasadena  1952  earthquake  record  was 
applied  to  the  bedrock  normalized  to  the  0.05g  level.  A  series  of  studies  were 
conducted  during  which  the  upper  layers  was  assumed  composed  of  sand 
corresponding  to  a  specific  Standard  Penetration  Test  blowcount.  The  blowcount 
was  converted  to  shear  modulus  and  the  associated  property  used  in  the  analysis. 
Figure  42  is  taken  from  DM7  and  is  based  on  work  by  Ohsaki  and  Iwasaki, 
Reference  28  from  1972.  This  data  is  old  and  limited  and  should  be  replaced. 
Reference  29  is  a  new  comparison  and  presents  Figure  43  which  shows  other 
relations.  For  the  purpose  of  this  analysis  the  most  representative  average  estimate 
was  thought  to  be  the  Ohta  relationship  for  sands.  Figure  44  shows  the  acceleration 
profiles  for  various  blow  counts  of  the  upper  70  foot  zone.  It  is  clear  that  the  looser 
the  zone  the  lower  the  shear  modulus  and  the  higher  the  amplification.  Figure  45 
illustrates  the  sharp  increase  in  amplification  when  the  blow  count  falls  below  20. 

Larger  Earthquakes  And  Treasure  Island 

A  question  occurs  as  to  how  Treasure  Island  will  respond  to  a  larger  level  of 
ground  shaking.  Will  there  be  a  corresponding  amplification  in  ground  motion?  To 
answer  this  the  analysis  described  above  was  repeated  for  an  earthquake  at  twice 
the  amplitude  of  the  Loma  Prieta  event.  Figure  46  shows  the  expected  site 
acceleration  with  depth.  Doubling  the  bedrock  motion  did  not  double  the  surface 
acceleration.  Nonlinear  behavior  resulting  in  a  softer  shear  modulus  combined  with 
increased  damping  reduced  the  amplification  from  about  3  to  about  2.  The  effect  of 
the  stiff  Bay  Muds  at  low  strain  produces  maximum  amplification  at  the  low  levels 
of  ground  shaking  associated  with  those  low  strains.  Thus  the  phenomenon  is 
primarily  associated  with  distant  earthquakes  having  low  levels  of  ground  motion 
which  excite  the  soil  deposits  at  low  levels  of  strain  to  produce  the  amplification. 
This  phenomenon  can  be  predicted  by  cyclic  shear  tests  to  establish  the  shear 
modulus  of  the  clays  under  various  levels  of  loading. 
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Figure  47  gives  guidance  based  on  Reference  31  for  estimating  the  effect  of 
Plasticity  Index  on  the  shear  stiffness  of  clays.  This  figure  is  suggested  for  use  on 
Navy  coastal  sites  where  high  plasticity  soils  like  Bay  Muds  can  be  found. 

SETTLEMENTS  OF  LIQUEFIED  SAND  DEPOSITS 

As  noted  above  settlements  occurred  on  Treasure  Island  and  caused  damage. 
Had  the  Loma  Prieta  earthquake  lasted  longer  it  would  have  induced  additional 
cycles  of  loading  in  the  soil  which  would  have  significantly  increased  the 
deformation.  As  a  saturated  cohesionless  soil  in  loose  condition  undergoes  cyclic 
shear  pore  water  pressures  begin  to  increase.  As  the  pore  pressures  increase,  the 
effective  confinement  is  reduced,  which  in  turn  causes  a  reduction  in  shear  modulus 
that  makes  it  easier  to  deform  the  soil  for  a  given  level  of  loading.  Unfortunately 
this  mechanism  feeds  upon  itself  such  that  as  the  soil  softens  with  reducing  shear 
modulus,  large  deformations  can  occur  with  the  onset  of  liquefaction.  Soils  which 
exhibit  restrained  deformations  under  one  level  of  shaking  can  exhibit  substantially 
greater  deformation  under  an  event  slightly  longer  in  duration  or  higher  in 
amplitude.  Navy  bases  situated  on  loose  saturated  soils  are  highly  vulnerable  to 
damage  unless  the  structures  have  been  designed  to  mitigate  this  effect.  Mitigation 
whether  in  the  form  of  ground  improvement  or  use  of  piles  is  very  expensive.  While 
we  can  predict  the  occurrence  of  liquefaction  on  level  ground  away  from  the 
structure  with  a  high  enough  degree  of  certainty,  prediction  of  liquefaction  beneath 
and  around  the  structure  is  substantially  more  complicated  and  not  adequately 
defined.  Prediction  of  the  deformation  state  (settlements  and  ground  flow)  is  even 
more  undefined.  This  is  a  major  area  of  concern  for  the  Navy  since  we  have 
numerous  existing  buildings  built  before  liquefaction  was  recognized  as  a  problem. 
New  construction  is  often  forced  to  include  the  most  conservative  foundation  from  a 
lack  of  understanding  of  the  liquefaction  deformation. 

Liquefaction  deformation  includes  two  basic  components:  settlements 
resulting  from  compaction  of  loose  saturated  deposits  with  an  outflow  of  pore  water 
under  shear  induced  pressures,  and  lateral  spreading  of  the  ground  from  flow  of 
reduced  strength  liquefied  deposits  associated  with  minor  slopes.  The  estimation  of 
deformation,  while  a  topic  for  research  by  many,  is  lacking  reliable  accurate 
procedures. 
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Settlement  Estimation 


The  state  of  the  art  for  settlement  prediction  of  liquefied  deposits  consists 
basically  of  empirical  data  relating  shear  strain  with  volumetric  strain.  In  Reference 
32,  Tokimatsu  and  Seed  present  data  shown  in  Figure  48  which  relates  the  volume 
strain  induced  by  a  magnitude  7.5  earthquake  as  a  function  of  effective  cyclic  shear 
strain.  They  suggest  that  other  magnitudes  may  be  scaled  by  use  of  data  in  Figure 
49.  The  authors  apply  this  data  to  an  estimation  of  settlement  shown  in  Figure  50. 
They  divide  the  soil  deposit  into  layers,  they  estimate  the  shear  modulus  based  on 
the  ratio  of  shear  at  a  strain  level  to  maximum  shear  strain,  and  then  use  relative 
density  or  blowcount  information  to  compute  volumetric  strain  using  Figure  48. 
The  volumetric  strain  is  adjusted  to  the  earthquake  magnitude  based  on  Figure  49 
and  then  corrected  for  multi-directional  effects.  Settlements  are  estimated  from  the 
volumetric  strain  of  each  layer  and  totaled. 

The  basic  procedure  outlined  by  Tokimatsu  and  Seed  can  be  substantially 
simplified  if  a  site  analysis  using  the  program  SHAKE  is  used.  Shake  computes  the 
effective  shear  strain  as  illustrated  in  Figure  51  for  the  Treasure  Island  site.  The 
calculations  can  be  improved  by  automation  in  the  form  of  a  spread  sheet.  A 
program  was  prepared  ,  LIQSS  which  performs  the  settlement  calculation  using  data 
from  Reference  32.  Figure  52  shows  the  spread  sheet  for  the  example  shown  in 
Figure  50.  The  computed  settlement  of  4.24  inches  for  a  magnitude  7.5  event  is 
reduced  by  0.8  for  a  magnitude  6.6  event  and  is  3.22.  The  value  of  3.22  is  slightly 
less  than  the  3.37  inches  computed  by  Reference  32.  The  difference  is  that  LIQSS 
uses  a  finer  interpolation  algorith  for  magnitude  correction. 

Data  from  the  SHAKE  analysis  of  Treasure  Island  was  used  with  blowcount 
data  of  the  site,  Figure  37  to  predict  the  settlements  at  location  UM10  on  Treasure 
Island.  The  results  are  shown  in  Figure  53  using  the  low  estimates  of  blowcount. 
Figure  54  using  average  estimates  of  the  blowcount  and  Figure  55  using  high 
estimates  of  the  blow  count.  Magnitude  corrected  settlements  range  from  4.8  inches 
to  12.7  inches  with  a  value  of  6  inches  for  average  values.  This  value  is  about  the 
same  as  the  6  inches  estimated  from  Figure  37. 

A  procedure  has  been  developed  using  the  data  from  Reference  32  using  the 
Tokimatsu  and  Seed  concept  but  modified  to  use  the  output  from  SHAKE  .  The 


procedure  gave  the  correct  order  of  magnitude  of  settlements  for  Treasure  Island 
when  average  blowcounts  were  used.  However  the  blowcount  data  is  imprecise  and 
the  range  of  possible  values  using  upper  and  lower  estimates  is  large.  For  this 
reason  NCEL  is  sponsoring  a  research  study  at  the  Waterways  Experiment  Station 
to  bring  together  a  group  of  experts  to  evaluate  the  available  data  and  attempt  to 
refine  the  process.  The  LIQSS  computer  program  is  an  interim  solution  to 
estimating  settlements  until  better  procedures  can  be  developed. 

CONCLUSION 

This  report  has  presented  a  detailed  study  of  the  response  of  the  Treasure  Island 
site  during  the  Loma  Prieta  earthquake.  It  shows  that  the  behavior  during  the 
earthquake  was  as  might  be  expected  given  the  loose  soil  conditions.  An  extensive 
discussion  of  settlements  and  liquefaction  was  presented.  There  was  an  extensive 
discussion  of  site  improvements  techniques  used  at  Treasure  Island  and  the 
resulting  damage.  The  site  was  analyzed  using  the  computer  program  SHAKE 
which  was  shown  to  give  an  accurate  assessment  of  site  response  and  soil 
amplification.  Choice  of  appropriate  material  properties  is  essential  for  an  accurate 
assessment,  although  the  analysis  was  shown  to  be  rather  insensitive  to  shear 
velocities.  SHAKE  should  be  used  as  a  preliminary  screening  tool  since  it  does  not 
consider  material  softening  with  development  of  pore  water  pressures;  nor,  does  it 
include  any  2-dimensional  effects  such  as  sloping  bedrock  conditions.  Caution 
should  be  used  in  using  code  type  spectra  for  soft  site  conditions.  The  effect  of 
initial  static  shear  stress  from  foundations  was  investigated  where  it  was  shown  that 
for  loose  soils  the  initial  shear  reduces  strength  rather  than  increases  it  as  might  be 
expected  for  dense  sands. 

Using  the  data  developed  by  Tokimatsu  and  Seed  an  automated  procedure 
was  developed  for  estimating  liquefaction  induced  settlements.  The  procedure  w  is 
applied  to  the  Treasure  Island  site  using  the  Loma  Prieta  earthquake  and  results 
agree  with  observed  settlements.  However  the  large  spread  in  blowcount  data  give 
a  large  range  in  estimated  settlements.  This  problem  is  under  study. 


Significant  findings  are; 


21 


•  Buildings  in  areas  which  settled  less  than  6  inches  generally  experienced  only 
minor  damage. 

•  For  loose  sands  the  initial  static  shear  caused  by  foundations  results  in  a  reduction 

of  strength. 

•  Site  improvement  at  least  during  moderate  events  was  very  effective. 
Alternative  procedures  appear  to  have  performed  equally  well. 

•  SHAKE  can  predict  site  response  and  soft  site  amplification  given  an  accurate 
assessment  of  the  subsurface  soil  properties. 

•  Treasure  Island  amplification  was  shown  to  be  dependent  upon  the  Bay  Mud 
properties. 

•  San  Francisco  Bay  Muds  resembles  the  Mexico  City  clays  in  being  stiffer  with 
straining  than  otherwise  expected.  This  is  related  to  its  Plasticity  Index 

•  Loose  sands  cause  site  amplification  especially  when  blow  counts  of  less  than  20 
occur  in  the  upper  layers. 

•  The  LIQSS  program  using  data  from  Tokimatsu  and  Seed  gives  reasonable 
estimates  of  settlement  when  average  blowcount  data  is  used. 

RECOMMENDATIONS 

The  following  are  recommended: 

•  That  Navy  review  soil  borings  for  bases  in  seismically  active  regions  to  identify 
regions  where  clay  deposits  are  prevalent. 

•  That  the  shear  modulus  properties  of  the  clay  under  various  levels  of  cyclic 
loading  be  obtained  from  soils  reports  if  available.  If  not  available  tests  should 
be  conducted  to  determine  whether  those  clays  exhibit  high  stiffness  as  found  in 
San  Francisco  and  Mexico  City.  Screening  should  include  evaluation  of  Plasticity 
Index 
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•  That  Figure  47  be  used  for  high  plasticity  clays  at  Naval  waterfront  sites. 

•  That  Navy  support  deep  instrumentation  arrays.  Such  an  array  will  yield 
information  on  actual  bedrock  motion  and  amplification.  NCEL  and  Western 
Division  have  jointly  supported  and  assisted  researchers  from  Brigham  Young 
University  and  the  University  of  New  Hampshire  perform  borings  and  field  tests 
at  Treasure  Island.  NSF  is  providing  funding  for  installation  of  a  deep 
instrument  array  on  Treasure  Island. 

•  That  Navy  pursue  NSF  funding  for  an  array  in  the  San  Diego  area. 

•  LIQSS  should  be  used  as  an  interim  approximate  means  for  estimating 
settlements.  Upper  and  lower  bound  blowcount  should  be  used  to  give  an 
estimation  of  the  range  of  possible  settlements. 
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Figure  1.  Location  of  Loma  Prieata  earthquake  and  aftershocks.  (Reference  2) 
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Figure  2.  Illustration  of  fault  movement  (Reference  2) 
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Figure  3.  Ratio  of  vertical  to  horizontal  acceleration,  both  components  considered. 
(Reference  1) 
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Figure  4.  Ratio  of  vertical  to  horizontal  acceleration,  average  of  horizontal 
components.  ( Reference  1) 
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Figure  5.  Map  of  affected  area.  (Reference  2) 
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Figure  6.  Map  of  active  faults  in  San  Francisco  Bay  area 
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Figure  7.  SanFrancisco  Bay  showing  Treasure  Island  and  soil  liquefaction. 
(Reference  4.) 
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Figure  9.  Treasure  Island  Site  plan.  (Reference  4) 


Figure  11.  Thickness  of  fill  and  native  shoals  materials.  (Reference  4) 
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Figure  13.  Treasure  Island  and  Yerba  Buena  showing  boring  locations. 
(Reference  5) 
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Figure  14.  Shear  wave  velocity  in  fill  and  shoal  sands.  (Reference  5) 
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Figure  16.  Normalized  shear  modules  and  damping  with  strain.  (Reference  5) 
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Figure  17  Treasure  Island  uncorrected  accelerogram.  (Reference  1) 
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Figure  18.  Yerba  Buena  uncorrected  accelerogram.  (Reference  1) 
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Figure  21.  Shaking  induced  compaction  settlement.  (Reference  4) 
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Figure  23.  Subsurface  profile,  Medical/dental  Clinic.  (Reference  13) 
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Figure  24.  Subsurface  profile,  Building  450  (Reference  13) 
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Figure  25.  Standard  Penetration  Test,  Building  450.  (Reference  13) 
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Figure  26.  Subsurface  profile.  Buildings  487, 488  and  489  (Reference  13) 
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Figure  27.  Subsurface  profile,  Building  453.  (Reference  13) 
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igure  28.  Spectral  accelerations  using  variety  of  assumptions.  (Reference  5) 
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SPECTRAL  ACCELERATION,  g  SPECTRAL  ACCELERATION. 


Figure  31.  Comparison  between  calculated  and  recorded  spectra  for  several  sites. 
(Reference  2) 
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(Reference  20) 
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Figure  34.  Folsom  gravel  Ka  (Reference  17) 
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Figure  35.  Variation  with  relative  density.  (Reference  20) 
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Figure  36  Types  of  contractive  deformation.  (Reference  20) 
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Figure  39  NCEL  computed  results  shown  with  results  from  Reference  5. 
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Figure  40.  Shear  modulus  as  function  of  strain. 
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Figure  41.  NCEL  computed  results  without  amplification  compared  with  Figure  39. 
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Figure  42.  Shear  modulus  as  function  of  blowcount. 
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Figure  43.  Shear  wave  velocity  as  function  of  blowcount.  (Reference  29) 
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Figure  48.  Prediction  of  volume  strain  from  cyclic  shear  strain.  (Reference  32) 
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Figure  50.  Example  study  from  Reference  32. 
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Figure  51.  SHAKE  analysis  results  of  Treasure  Island  site  Loma  Prieta  earthquake. 


=|  READY  k| 


D28  6.6 


A 

B 

c 

D 

E 

F 

1 

2 

3 

4 

Tokimatsu  Seed  Liquefaction 

Settlement 

Layer 

Thickness 

BlowCount 

Eff  Shear 

Vol  Strain 

Settlement 

5 

- ft  — 

1  —  — — — 

—— —  %  — — — — 

- In - 

6 

1 

5.00 

9.00 

0.0500000 

0.1460495 

0.18 

7 

2 

5.00 

9.00 

0.0800000 

0.2336792 

0.28 

8 

3 

10.00 

9.00 

0.1200000 

0.3505188 

0.84 

9 

4 

10.00 

9.00 

0.1400000 

0.4089387 

0.98 

10 

5 

10.00 

9.00 

0.1500000 

0.4381486 

1.05 

11 

6 

10.00 

9.00 

0.1300000 

0.3797288 

0.91 

12 

7 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

13 

8 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

14 

9 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

15 

10 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

16 

11 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

17 

12 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

18 

13 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

19 

14 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

20 

15 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

21 

16 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

22 

17 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

23 

18 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

24 

19 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

25 

20 

0.00 

0.00 

0.0000000 

0.0000000 

0.00 

26 

27 

28  Enter  earthquake  magnitude  6.60 

29 

30 

31 

32  Total  settlement  in  inches  3.22 

33 

34  Enter  data  in  columns  B  C  and  D  for  each  sand  layer 

35  Column  B  Layer  thickness  in  Feet 

36  Column  C  Blowcounts  N  sub  I 

37  Column  D  Effective  Shear  Strain  in  percent 

38 

39 


What  would  like  to  do?  P(rint),  S(ave  File),  G(et  File),  Q(uit) 


Figure  52.  LIQSS  results  for  example  in  Figure  50. 
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Figure  53.  LIQSS  results  for  Treasure  Island  using  lower  bound  blowcounts. 
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Figure  54  LIQSS  results  for  Treasure  Island  using  average  blowcounts. 
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Figure  55.  LIQSS  results  for  Treasure  Island  using  upper  bound  blowcounts. 
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Appendix 


Computer  Program  LIQSS 


INTRODUCTION 

The  computer  program  LIQSS  was  developed  as  an  interim  solution  for 
computing  settlements  from  liquefaction  of  cohesionless  soils.  It  is  based  on  data 
from  Reference  32 .  The  procedure  requires  a  SHAKE  analysis  be  perfromed  for 
the  site. 

SYSTEM  REQUIREMENTS 

The  program  LIQSS  will  operate  on  any  PC  computer  with  a  minimum  512  k 
memory  and  DOS  3.2  or  higher.  Printing  of  results  is  accomplished  by  a  parallel 
printer  connected  to  LPT1  or  a  serial  printer  with  LPT1  redirected. 

INPUT  REQUIREMENTS 

A  SHAKE  analysis  is  required  for  use  of  LIQSS  program.  SHAKE  requires 
the  soil  deposit  be  divided  into  layers,  and  the  design  earthquake  time  history  record 
be  specified,  and  shear  and  damping  properties  of  the  soil  layers  be  specified.  The 
user  is  expected  to  be  familiar  with  this  procedure. 

In  addition  to  the  SHAKE  parameters,  blowcounts  data  in  the  sand  layers  is 
required.  Place  the  disk  in  drive  A  and  type  LIQSS  to  begin  the  program.  A  spread 
sheet  will  be  presented.  Enter  data  for  each  sand  layer.  Do  not  include 
nonliquefiable  layers.  Enter  the  effective  cyclic  shear  strain  from  the  SHAKE 
analysis  results.  As  each  layer  is  completed  the  settlement  is  computed  and  the  total 
displayed.  To  save  or  print  results  press  the  slash  key ,  /;  then  enter  P  to  print,  S  to 
save  or  G  to  get  a  rile  from  disk.  Note  the  spread  sheet  extends  from  colums  A 
through  F  and  rows  1  through  37.  Depending  on  the  type  of  monitor,  not  all  of  the 
screen  may  be  visible.  Use  the  arrow  keys  to  move  about  the  spread  sheet. 

The  following  definitions  apply: 


A-l 


Thickness  Thickness  of  a  layer  in  feet.  Use  same  layers  used  in  SHAKE. 

Blowcounts  Use  Nj  corrected  blowcount  for  layer. 

Effective  Use  SHAKE  data  expressed  as  percent. 

Shear  Strain 

Magnitude  Richter  magnitude  of  earthquake 

Computed  results  are: 

Volume  Strain  Volume  Strain  in  percent 

Settlement  Settlement  of  layer  in  inches 

Total  Settlement  Total  settlement  for  all  liquefiable  layers  in  inches. 


A- 2 


DISTRIBUTION  QUESTIONNAIRE 

The  Naval  Civil  Engineering  Laboratory  is  revising  its  primary  distribution  lists. 


SUBJECT  CATEGORIES 

1  SHORE  FACILITIES 

1A  Construction  methods  and  materials  (including  corrosion 
control,  coatings) 

1 B  Waterfront  structures  (maintenance/deterioration  control) 

1C  Utilities  (including  power  conditioning) 

1 0  Explosives  safety 
IE  Aviation  Engineering  Test  Facilities 
1 F  Fire  prevention  and  control 
1G  Antenna  technology 

1 H  Structural  analysis  and  design  (including  numerical  and 
computer  techniques) 

1J  Protective  construction  (including  hardened  shelters,  shock 
and  vibration  studies) 

1 K  Soil/rock  mechanics 
1 L  Airfields  and  pavements 
1 M  Physical  security 

2  ADVANCED  BASE  AND  AMPHIBIOUS  FACILITIES 

2A  Base  facilities  (including  shelters,  power  generation,  water 
supplies) 

2B  Expedient  roads/airfields/bridges 
2C  Over-the-beach  operations  (including  breakwaters,  wave 
forces) 

2D  POL  storage,  transfer,  and  distribution 
2E  Polar  engineering 

3  ENERGY/POWER  GENERATION 

3A  Thermal  conservation  (thermal  engineering  of  buildings, 
HVAC  systems,  energy  loss  measurement,  power 
generation) 

3B  Controls  and  electrical  conservation  (electrical  systems, 
energy  monitoring  and  control  systems) 

3C  Fuel  flexibility  (liquid  fuels,  coal  utilization,  energy  from  solid 
waste) 


3D  Alternate  energy  source  (geothermal  power,  photovoltaic 
power  systems,  solar  systems,  wind  systems,  energy 
storage  systems) 

3E  Site  data  and  systems  integration  (energy  resource  data, 
integrating  energy  systems) 

3F  EMCS  design 

4  ENVIRONMENTAL  PROTECTION 

4A  Solid  waste  management 

4B  Hazardous/toxic  materials  management 

4C  Waterwaste  management  and  sanitary  engineering 

4D  Oil  pollution  removal  and  recovery 

4E  Air  pollution 

4F  Noise  abatement 

5  OCEAN  ENGINEERING 

5A  Seafloor  soils  and  foundations 

5B  Seafloor  construction  systems  and  operations  (including 
diver  and  manipulator  tools) 

5C  Undersea  structures  and  materials 
5D  Anchors  and  moorings 

5E  Undersea  power  systems,  electromechanical  cables,  and 
connectors 

5F  Pressure  vessel  facilities 
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5H  Ocean-based  concrete  structures 

5J  Hyperbaric  chambers 

5K  Undersea  cable  dynamics 
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